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ABSTRACT 


ERROR  ANALYSIS  OF  CASH  METHOD:  COMPUTER-AVERAGING 


OF  NOISY  PERIODIC  SIGNALS 


Samuel  A.  Elder 

Physics  Deiiartment 
United  States  Naval  Academy 
Annapolis,  Maryland  21402 


Computer-assisted  synchronous  hot  wire,  or  "CASH,"  method 
was  developed  to  permit  computer  averaging  of  hot  wire  data  when 
sampling  oscillating  velocity  profiles.  Using  emitted  sound 
pressure  as  phase  reference,  computer  separates  periodic  component 
of  hot  wire  signal  from  turbulent  background  by  averaging  over 
large  number  of  acoustic  cycles.  So  far  technique  has  been  utilized 
in  two  investigations:  (1)  interface  oscillations  of  flow-excited 

cavities,  and  (2)  growth  of  instability  waves  in  laminar  jets. 
Accuracy  depends  on  the  ability  to  match  closely  the  reference  source 
frequency.  Any  mismatch  produces  systematic  error.  For  given 
amount  of  mismatch,  however,  error  is  minimized  by  dividing  sample 
into  subsets  and  performing  the  averaging  at  two  levels.  In  the 
cavity  resonator  problem,  random  deviations  from  periodicity  in 
the  reference  itself  make  frequency-matching  difficult.  Measured 
standard  deviations  show  a predicted  periodicity,  with  greatest 
deviation  occurring  90°  after  signal  maximum.  By  counting  beats 
against  a standard  oscillator,  frequency  error  can  be  estimated  at 
data  acquisition  time.  With  this  information,  sync  error  in  average 
signal  can  be  completely  removed  by  post-run  data  reduction. 


Work  sponsored  by  Naval  Ship  System  Command,  GHR  program,  administered 
by  Naval  Ship  Research  and  Development  Center. 
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ERROR  ANALYSIS  OF  CASH  METHOD:  COMPUTER- AVERAGING 

OF  NOISE  PERIODIC  SIGNALS 
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Samuel  A.  Elder 

United  States  Naval  Academy 
Annapolis,  Maryland  21402 

Slide  I. 

The  computer-assisted  synchronous  hot  wire,  or  "CASH"  method  was 
developed  to  permit  computer  averaging  of  hot  wire  data  when  sampling 
oscillating  velocity  profiles.  In  the  figure  a cavity  is  shown  being 
excited  by  fluid  flow.  A hot  wire  probe  scans  the  field  in  the  mouth 
opening,  recording  velocity  data  in  computer  core.  Simultaneously 
the  signal  from  the  sound  pressure  microphone,  inside  the  cavity,  is 
also  recorded,  to  give  a phase  reference. 

Later,  the  computer  data  array  is  averaged  to  remove  turbulence 
noise,  and  re-arranged  to  give  the  time  sequence  of  vertical  profile 
oscillations  in  the  shear  layer. 

Slide  2. 

The  idea  is  to  separate  the  periodic  component  of  the  hot  wire 
signal  from  background  turbulence.  Turbulence  noise,  being  random,  is 
easily  removed  by  averaging,  but  to  preserve  the  phase  relation  in  the 
pressure  signal,  the  averaging  must  be  done  synchronously.  Typical 
initial  noise  level  is  100%  of  signal,  i.e.,  S/N~l.  After  10  periods 
of  averaging  N/S  level  is  reduced  by  factor  3,  after  100,  by  factor  10, 
etc.,  since  noise  drops  as  1/i/M. 

Length  of  sample  run  was  not  extended  beyond  100,  due  to  core  size 


limitation . 
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Slide  3. 

In  addition  to  vertical  velocity  profiles,  it  is  possible  to  form 
horizontal  profiles  of  the  interface  displacement  oscillation,  by 
a further  reduction  of  the  data. 

Here  is  a sequence  of  interface  wave  motions  for  a cavity  oscilla- 
ting at  283  llz,  taken  at  0°  (i.e.,  in-phase  with  pressure),  90°,  180°, 
and  270°  in  the  cycle.  Flow  is  to  the  right  at  a free  stream  speed  of 
27  m/s,  with  leading  edge  at  left,  and  trailing  edge  at  right.  Circles 
represent  planes  at  which  measurement  was  performed. 

Slide  4. 

In  order  to  obtain  a meaningful  average,  the  periodic  component 
of  the  velocity  signal  is  sampled  at  exactly  20  times  the  sound 
frequency  for  100  periods  of  the  acoustic  signal.  In  principle  this 
should  give  a faithful  reproduction  of  the  periodic  part  of  the 
velocity  signal,  out  to  the  10th  harmonic. 

To  determine  the  reliability  of  the  method,  the  computer  was  pro- 
grammed to  calculate  the  normalized  standard  deviation  of  the  velocity 
signal  during  the  averaging  process.  The  measured  standard  deviation 
when  plotted  as  a function  of  time,  was  found  to  have  a large  periodic 
component,  suggesting  some  systemmatic  error  process  involved  in 
the  method.  A likely  source  was  suspected  to  be  the  inability  of  the 
experimenter  to  obtain  perfect  synchronism  between  the  sound  and 
sampling  frequencies. 

Slide  5. 

Suppose  that  there  is  an  error  Af  in  setting  the  sampling  frequency 
f,  then  there  will  be  an  imperfect  phase  closure,  AO^  in  each  cycle. 
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causing  the  apparent  period  to  be  given  by  T'=T(l+Af/f)  as  shown. 

The  measured  signal  level  at  phase  0 in  the  ith  cycle  will  be 

A.  = A cosfO-({)  +A0 . ) + e 

1 p ^ ^1  i-" 

where  A^  = "True"  amplitude  of  periodic  signal 

d)  = initial  offset 
1 

e = random  noise  component 
Slide  6. 

Computing  the  average  over  M-periods,  we  see  that,  for  small  total 
closure  error  A0|^j,  the  effect  of  sync  error  cancels  out,  leaving  only 
the  "true"  signal  after  averaging.  (The  random  error,  c,  by  hypothesis, 
is  taken  to  average  out  to  zero.) 

Slide  7. 

From  the  estimated  mean,  the  variance  can  be  formed,  as  shown. 
Again,  for  small  the  expression  may  be  simplified  somewhat. 

The  residual  variance  has  both  periodic  and  random  components. 

Slide  8. 

The  Standard  Deviation,  normalized  by  the  amplitude  of  the  first 
Fourier  Component  of  the  signal,  is  evaluated  by  summing  over  the 
geometric  series  in  A0^,  The  periodic  component  of  o is  found  to  vaVy 
with  M.  (For  large  M,  the  variation  approaches  simple  proportionality.) 

Slide  9. 

Here  is  what  ..ic  com]nitcr-averaged  data  and  the  NSI)  look  like  when 
plotted  over  one  period.  Data  are  represented  by  circles,  while 


solid  lines  represent  theoretical  predictions.  The  frequency  error, 

Af,  estimated  by  counting  beats  between  the  sound  and  the  standard 
oscillator  signal  was  found  to  be  about  1.2  Hz. 

Note  the  distinctly  periodic  nature  of  tlie  error.  Fitting  the 
deviation  curve  to  the  theoretical  model  was  done  by  taking  the  value  at 
0 = to  be  the  random  noise  component. 

Despite  the  large  deviation  near  90°  and  270°,  the  fit  to  the 
sampled  sound  pressure  was  very  close,  as  predicted.  In  this  run 
Sound  Pressure  was  analyzed  rather  than  velocity,  because  its  nearly 
sinusoidal  waveform  best  fitted  the  theoretical  model. 

Slide  10. 

Before  dealing  with  the  velocity  data,  the  more  general  case  must  be 
treated.  This  is  done  by  expressing  the  signal  as  a Discrete  Fourier 
Series  of  N terms,  where  N = '2  the  number  of  points  sampled  per  period. 
Also  the  requirement  of  small  closure  error  is  dropped.  IViien  the 
average  is  taken,  e again  vanishes  by  hypothesis,  while  sinA6^  is 
found  to  vanish  due  to  symmetry.  The  cosA0^  term  turns  out  to  be  the 
familiar  single-slit-diffraction  function,  sin  x/x.  Since  this  term 
remains  as  a "weighting  function"  in  the  average-value  series,  we  see 
that,  in  the  general  case,  the  effect  of  sync  error  does  not  automatically 
cancel  out! 

Slide  11. 

Here  the  error-weighting  function  is  plotted  for  several  values 

of  M. 

For  a given  value  of  a,  (and  thus  of  Af/f) , least  error  is 
produced  in  the  answer  by  averaging  all  the  data  in  pairs  of  cycles, 
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i.e.,  taking  M=2.  Unfortunately  this  requires  maximum  computer  processing 
time.  For  practical  puriwscs,  a value  of  M=10  was  found  to  be  most 
satisfactory. 


Slide  12. 

Here  is  shown  the  computation  of  normalized  standard  deviation 
for  the  general  case.  Despite  cancellations,  the  expression  for 
o/A^^j  remains  complicated.  Dependence  on  Af/f  can  be  illustrated  by 
performing  the  average  of  sin^AB.  numerically  for  several  values  of  M. 


Slide  13. 

The  calculations  are  shown  here. 
Af/f,  the  error  rises  more  quickly  as 
the  maximum  value  of  sin^AB.  actually 


It  may  be  seen  that,  for  small 
M gets  larger.  However, 
occurs  for  M = 2. 


Slide  14. 

From  the  variation  of  cosABj  with  M and  a,  one  can  see  why  a value 
of  ten  is  a good  choice  for  M.  An  arrow  shows  the  abscissa  of  the 
highest  harmonic  for  the  data  analyzed  (i.e.,  the  10th).  Sync  error 
for  this  term  is  no  greater  than  15%  at  M = 10,  while  for  lower 
harmonics  the  error  is  considerably  less.  Since  100  periods  must  be 
averaged  to  properly  reduce  the  background  noise,  a choice  of  M = 10 
means  that  the  data  must  be  grouped  into  ten  subsets  of  ten  periods  eacli. 

If  M = 2 had  been  chosen,  the  error  would  have  been  less,  but  the 
computation  time  would  have  been  increased  by  a factor  of  nearly  five  times, 
since  50  subsets  would  be  formed.  On  the  other  hand,  if  M = 100  were 
chosen,  performing  the  average  in  a single  sweep,  gross  distortion  of 
the  signal  would  result.  Since  the  exact  form  of  the  weighting  factor  is 
known,  an  interesting  possibility  arises  that  the  sync  error  could  be 


completely  removed  by  a post-run  calculation,  provided  a reliable 
estimate  of  sync  frequency  error  is  obtained  at  data  acquisition  time. 
(The  beat-counting  method  used  previously  would  be  suitable  for  this.) 
However,  because  of  the  additional  computational  time  that  this  would 
introduce,  it  is  more  convenient  to  choose  an  operating  region  for  which 
the  sync  error  can  be  ignored  relative  to  ether  experimental  errors. 


6 


CAVITY  RESONATOR 


APPARENT  PERIOD: 


T'  = T-C  1 + Af/f} 


PHASE  ERROR  AT 

ith  CYCLE:  «•  E tt  i Af/f 

INSTANTANEOUS  VALUE:  A.  = A cosI0  -6  ♦ A6  } + e 

1 p i-^ 

UHERE: 

f=  ESTIMATED  SIGNAL  FREQUENCY 
Af=  ERROR  IN  f 

Ap=  "TRUE"  AMPLITUDE  OF  PERIODIC  SIGNAL 
INITIAL  OFFSET 
e=  RANDOM  NOISE  COMPONENT 


Fig  S 


EXPANDING-. 

A = A cos-C9-(().>  cos  A0  - A sirvC0-((),>  sin  A0„  7 

P I n p 1 n 


STANDARD  DEVIATION 


= |A.  - A(  ^ 

SUBSTITUTING-, 

^ S S E ~ 

o = Ap  sin  -Ce-^jIAei,^  - EeAO^  A^  sin{0-(t)j}  + e 

OR 

^ 2 . E E 2 

a - sin'^-Ce-(j),>  A0„  ■*■  a 

p 1 N n 

UHERE  = Std.  Dev-  of  random  component 


Fig  7 
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Fig  10 
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NORMALIZED  STANDARD  DEVIATION 


Fig  m 
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oscillating  velocity  profiles.  Using  emitted  sound  pressure  as  phase 
reference,  computer  separates  periodic  component  of  hot  wire  signal 
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(1)  interface  oscillations  of  flow-excited  cavities,  and  (2)  growth  of 
instability  waves  in  laminar  jets.  Accuracy  depends  on  the  ability  to 
match  closely  the  reference  source  frequency.  Any  mismatch  produces 
systematic  error.  For  given  amount  of  mismatch,  however,  error  is 
minimized  by  dividing  sample  into  subsets  and  performing  the  averaging 
at  two  levels.  In  the  cavity  resonator  problem,  random  deviations  from 
periodicity  in  the  reference  Itself  make  frequency-matching  difficult. 
Measured  standard  deviat.ions  show  a predicted  periodicity,  with  greatest 
deviation  occurring  90^^  after  signal  maximum.  By  counting  beats  against 
a standard  oscillator,  frequency  error  can  be  estimated  at  data 
acquisition  time.  With  this  information,  sync  error  in  average  signal 
can  be  completely  removed  by  post-run  data  reduction. 
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